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Fonctions de Bessel 1 et l'équation du temps de Kepler 2 


Corrigé par M.TARQI 

Parie I. Une fonction définie par une intégrale 


1.1 Si À > 0, la fonction 
Si À < 0, on a 


Q —ix cos 6 2À 


sin"’ y est continue sur [0, 7 r], donc elle est intégrable sur ]0, 7 r[. 

2A a 


g— ix cos 0 sin 2A ( 


= 5111 " 


Donc il suffit de montrer que la fonction 9 \ — > sin 2A 9 est intégrable sur ]0, 7r[, en effet, tout fonction 
absolument intégrable est intégrable. La fonction 9 \ — > sin 2A 9 est continue sur ] 0, tt[, au voisinage 

de 0, sin 2A 9 


et au voisinage de -k 


(tt - 9)~ 2X sin 2A 9 = (tt - 9)~ 2X S ïn 2X (n - 9) = ( ^ ~ 

\ (.v-n) / 


2A 


1, 


et comme 0 < — 2À < 1, la fonction ( 
1.2 D'après la relation de Chasles, on a 


sin 


2A 


9 est intégrable sur ]0, 7r[. 


/a(x) = [ e 

Jo 


I 2 

— ixcosO gj n 2A Q d Q / ixcosti „-„2A 


j; 


sin 2A 9d9+ / e~ ix cos 6 sin 2X 9 d9 


Le changement de variable t = n — 9 montre que 


0 -ix cos e sin 2A Q de= _ / e ix cos t sin 2A ^ 


D'où 


f\(x) = 


ri\ 

J 0 


e- ixcos6 sm 2X 9d9 


—ix cos 6 _• 2A 


£■ 


sin ZA 9d9- I e ix cos t sin 2A tdt 

7T 
2 

,ix cos t\ _• 2 A . 


2Re(e“ cost )sin 2A tdt 


= 2 


fl 

/ cos ( x cos t ) sin 2 A tdt 

J o 


1. Bessel, Friedrich Wilhelm (1784-1846), astronome et mathématicien allemand, connu principalement pour avoir effectué 
les premières mesures précises de la distance d'une étoile et pour être le fondateur de l'école allemande d'astronomie d'ob- 
servation. Il introduit, dans la résolution des problèmes de mécanique céleste faisant intervenir la théorie des perturbations, 
les fonctions mathématiques dites de Bessel, solutions d'équations différentielles particulières. Ces fonctions jouent un rôle 
important dans l'analyse de la répartition et de la conduction de la chaleur ou de l'électricité à travers un cylindre. 

2. Johannes Kepler (1571-1630), est un astronome allemand célèbre pour avoir étudié l'hypothèse héliocentrique (la Terre 
tourne autour du Soleil) de Nicolas Copernic, et surtout pour avoir découvert que les planètes ne tournent pas en cercle parfait 
autour du Soleil mais en suivant des ellipses. 
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1.3 Posons / A (x) = / g\(x, 9)d9 avec g\(x, 9) = e~ zx cos 9 sin 2A 9. Montrons par récurrence sur p G IN, 

J o 

que fx est de classe Y? p sur IR. 

Pour p = 0, la fonction x i — > e~ lx cos 9 sin 2A 9 est continue sur IR pour tout 9 g]0, tt[, de plus 
\g\(x, 9)\ < 9 sin 2A 9 et 9 \ — > sin 2A 9 est intégrable, donc d'après le théorème du cours x i — > f\(x) 
est continue sur IR donc de classe Y A) sur IR. Supposons maintenant / est de classe Y' 7 ’ sur IR et que 

Vx € IR, f {p \x) = ( -i) p f* e ~ ixcose cos p 9 sin 2A 9d9. 

J o 

La fonction (x,6) i — > (—i) p e~ lx cos e cos p 6 s'm 2X 9 est continue sur 1 Rx]0,7t[ et admet une dérivée 
partielle continue par rapport à x, 

(x,9) ^ {-i) p+1 e- ixcosd cos p+1 9sm 2X 9 

qui est dominée par la fonction intégrable : 

(px : 9^ | sin 2A 9\. 

Donc, toujours d'après le théorème de cours, est Y’ 1 sur IR, c'est à dire fx est c é ,p+1 sur IR. 

Ainsi fx est de classe sur IR et : 

Vx € IR, fx’\x) = (—i) p [ e~ lx cos 6 cos p 9 s'm 2X 9 d9. 

J o 


1.4 Une équation différentielle 

1.4.1 Pour tout x G IR, on a : 

/a+i(x) - f x (x) = r e~ ixcos6 9( sin 2A+2 9 - sin 2X )d9 = - T e~ ix cos ' 9 cos 2 9 sin 2A 9d9. 

J o J 0 

Donc f" = fx+i - fx. 

1.4.2 Soit xGIR. 


-x/ a+ i(x) = i n e~ ix cose y sin 2X+1 9d9 

J o 


= % 


—ix cos 9 „:„2A+1 


sin + 9 


/*71 

J 0 


—ix cos 0 


d 


— (sin 2X+1 9)d9 


= i 


i{ 2A + 1) f 

J o 


J o J 0 d9 

r j—ix cos 0 sin 2A 9 cos Q d g = (2A + 1 )fx(x). 


D'où l'égalité demandée. 

1.4.3 Pour tout x G IR, on a : 

xf'Rx) + (2À + 1 )fx(x) + xfx(x) = x(/a+ i(x) - fx(x)) + (2A + l)/i(x) + x/ A (x) = 0 

Donc fx est solution de l'équation différentielle (F A ) sur IR. 

Parie II. Une formule d'Euler 


2.1 
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2.1.1 A l'aide d'une intégration par parties, on obtient : 

1 


r+oo 

r (p)= / t p ~ 1 e~ t dt = 

J 0 


-t p e 


—t 


+00 -j /*+oo 1 

+ - / t p e- l dt = -r(p+l). 


LP Jo PJO 


P 


r-\-oo 

2.1.2 On a T(l) = / e t dt = 1, donc la formule précédente entraîne que, pour tout n € IN*, 
2o 

T(n) = ni 


f+OO 


2.1.3 La fonction ti — > t p e 1 étant positive sur ]0, +oo[ donc T(p) = / t p e > 0. Supposons 

J o 

qu'il existe p > 0 tel que T (p) = 0, alors on aura pour tout a > 0, / t p ~ l e~ t = 0 et la continuité 

J o 

de t 1 — > t p ~ 1 e~ t implique = 0 pour tout t > 0 ce qui est absurde. 

2.1.4 II suffit d'utiliser le changement de variable t = s 2 . 


2.2 


2.2.1 La fonction de deux variables : (x, y) i — > x 2p 1 y 2 ' 1 1 e x y est une fonction à variables sépa- 
rables, donc par le théorème de Fubini, on peut écrire : 


[0,a,] ‘ 


x 2p-l y 2 q -l e -x 2 -y 2 


dxdy = f x 2p 1 e x2 dx f y 2q 1 e y2 dy = I p (a)I a (q ) 

J o Jo 


2.2.2 Utilisons les coordonnées polaires :(x, y) = (r cos 9 , r sin 9) . Quand ( x , y) décrit D(a), r décrit 


TT 


l'intervale [0, a] et 9 décrit l'intervalle ^0, —J . La formule de changement de variables s'écrit 
alors : 


[0,o] ■ 


2.2.3 II est clair que T(p) = lim I p (a), donc 


x 2p-i y 2q-i e -x 2 -y 2 dxdy = / / 2 cos 2p - 1 9sm 2q - 1 9r 2( P + ^- 1 e- r2 drd9 


0 J 0 


r 2(p+ ? )-l e -r 2 dr / 2 cos 2p-l Q gin 2g-l m 


Ip+ q (a ) / cos 2p 1 #sin 2,? 1 9d9. 

Jo 


lim / / x 2p l y 2q 1 e x2 yl dxdy = r(p)r(g). 
a ^°°J J [o , o] 2 

D'autre part, vues les inclusions entre les domaines d'intégrations D(a), [0, a] 2 ,D(aV 2) et vu 
le fait que V(x, y) G 1R 2 +, x 2p ~ 1 y 2q ~ 1 e~ x2 ~ y2 > 0, on a : 


x 2 p l y 2t ? 1 e x V dxdy < I p (a)I q (a ) < / / x -1 ' ‘y 1 *e 

D(a) J J D(ay/2) 


2p-l„2q-l P -x 2 -y 2 dxdy. 


Or 


lim / / x 2p 1 y 2lJ x e æ2 y2 dxdy = lim 

a^oo / / D(o) a^ooj J D(a y2) 


x 2p - l y 2q - l e- x -V dxdy 


= r (p + p) cos 2p 1 0sin 2q 1 
do 
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alors par comparaison des limites il vient : 

H 

T(p)T(q)=T(p + q) / 

J 0 

2.3 Posons u = \[t, alors : 


cos 2p_1 6 sin 29-1 QdO. 


a-r 


e f = Vvr 


D'autre part, la relation précédente dorme 


r - r(n) = r n + - / 2n _ 


L 


avec I p = sin p xdx. Pour tout p > 2, on écrit à l'aide d'une intégration par parties. 


Ip — 


f: 


= [- 


sin p 1 x sin xdx 

i n - ( ■ 

®]o +(?-!) / 
Jo 


cos x sm 


sin p 1 x cos 2 xdx 


= (p — 1) / sin p 1 x(l — sin 2 x)dx = (p — l)(/ p _2 — i p ) 
iO 

p — 1 7T 

d'où I p = Ip— 2 - En tenant compte de Iq = — et I i = 1, on a alors 

p 2 


lin - 1 — 


(2n — 2)(2n — 4). ..4. 2 2 n " 1 ((n-l)!) 

(2n — l)(2n — 3). ..3 ~~ (2 n - 1)! 


Ainsi 


F f n 4- — ^ = ^ 2n ~ X ) ! = ( 2n V-V^ 

V 2 J 2 n_1 (n — 1)! 2 2n n! 

Parie III. Fonction de Bessel 


3.1 On a, pour tout x e]0, +oo[, /i (x) = — 

2 IX 


Jl(x ) = 


-iarcoael 7 " 2 sin X „ 

e “ costf = , d où: 

J 0 X 


1 


r(è)r(i) V 2 2 


nfà(x) = \ ~ 


2 sinx 


vr Wx 


Donc Ji (x) 


^ y/x au voisinage de 0 + (lim = 1). 

TT x—>0 X 


3.2 En tenant compte de la relation — x/a+i(x) = (2À + l)f' x (x), on a, pour tout x €]0, +oo[ : 


A 


A 


1 


x Jx Jx xr(i)r(A + \) ( 2 ) fx ^ 2r(i)r(A + i) ( 2 ) fx ^ 


+ 


1 


X\ x x/a+i(x) 


r(i)r(A + i) \2) (2A + 1) 

1 /x\ a x/a+i(x) 


r(i)r(A + i) \2) (2A + 1) 

1 /x\ A +1 /a+i(x) 


r(è)E(A + i) V 2 J (A + i) 


1 


A+l 


P(è)E(A + §) V2 

= Ja+i(t)- 


/a+i(t), car r I A + - ) — I A + - A + 


1 



1 
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ld_ 

x dx 


3.3 Montrons par récurrence sur n G IN*, que pour tout / G ^°°(]0, +oo[, ]R) et tout x > 0 : 

[H n ~\f)Ux) = xH n {f)(x). 

La propriété est vraie pour n = 1 , puisque [H°(f)]' (x) = f\x) = X 
f G ^°°(]0, +oo[, IR), H(f) G <r°°(]0, +oo[, IR) et donc 

Vx > 0, [H n {f)]\x) = [H n - 1 H{f)}'{x) = xH n (H(f))(x) = xH n+ \f){x ) 

Ainsi, par le principe de récurrence, Mn G IN* 

V/ G tf°°(]0,+oc[,]R), Vx > 0, [H n -\f))'{x) = xH n if)(x). 


(/)(x) = xH(f){x). Si 


La propriété J n+ i(x) = (— l) n \ — x 


2 .n+è i\ d 


x dx 


smx 


x 


est vraie pour n = 0, d'après la ques- 


tion 3.1. Supposons qu'elle vraie à l'ordre n et montrons la pour n + 1. On a, grâce de la relation 

[H n - 1 (f)]'(x) = xH n (f)(x): 


J n + Î = 


n + \ 


x 


2 J n+ Ux)~ J' i(x) 


n+ ï ( - 1) n/l x n + i(lA. 

V n \ x dx 


1 d 




X 


-(-1 ) n \l-x n+ ïx 


1 d \ 


n+l 


x dx) 


smx 

x 

sinx 

x 

sinx 


= (-i) 


n+l 


,n+§ 


1 d \ 


n+l 


U ~xTx 


f smx 

V x 


- (=S 

_(_ irv /f(n + l)x”-> 


D'où l'égalité demandée. 

3.4 Pour tout x G]0, +oo[, on a, en notant A = r( + r(A + -4): 


smx 


1 d_ 

x dx 


smx 


A (x 2 J\(x) + xJ' x (x) + (x 2 — A 2 ) J A (x)) = x a A(A 4 ^ / A (x) + 2x 2 ^(|) / a (x) 


+x 


/"(*) + * 


A— 1 


f\(x) + X 


X 


/a (s) 


+(x 2 - A 2 )/ a (x) 


[A(A - 1)/ A (x) + 2 xA/ a (x) 

+x 2 +'(x) + A/ a (x) + x/ A (x) + (x 2 - A 2 )/ a (x)] 
= (f) [^ 2 /a (^) + (2A + 1)x/;(x)x 2 / a (x)] 

= x [x/ A (x) + (2A + 1)/ A (x)x/ A (x)] , 

et comme / A est solution de l'équation différentielle F\, alors 

Vx > 0, x 2 J a (x) + xJ A (x) + (x 2 — A 2 )J a (x) = 0, 
donc J A est solution de l'equation différentielle L A sur ]0, +oo[. 
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±2 n 


3.5 On a tout réel t, cos t = E(- 1 >”rânï- d ' oi,: 


n = 0 


oo / a \ 2)2 

Vx > O,V0 G],7r[, cos(xcos$) sin 2A 9 = ^(— l) n sin 2A — . 


n=0 


Posons u n (6) = (— l) n sin 2A 0 


n • 2A ,(æcos0) 2n 


(2n) 


pour tout n G IN. Les sont intégrables sur 


7 r 

°’2J 


puisque 


2n 

X ■ • 2A, 


|u n (0)| < t — — | sin 2A 0\ et / | sin 2A 0|d0 existe ( la question 1.1 ). La dernière inégalité montre que 

(277.) ! J o 

C'h 

la série YJ / \u n (0)\d0 est convergente. 

nelN Jo 

D'où, par le théorème du cours et en tenant compte de la relation (1) de la question 2.2.3 : 


f\(x) = 2 / cos(zcos0)sin 2A 0d0 

J o 


2 £(-!)' 


J2n 


n = 0 
oo 


(2n)! 


cos 2n 6 sin 2A 6d6 


B-i)^ r( - +è)r(A+è) <*> 


n= 0 


(2n)! r(n + À + l) 


et par conséquent, en utilisant l'égalité T ( n + ^ j = - ^ : 

Zj J Zj Tll 


J\{x) = 


r( i)r( n+ i) V 2 


r(è)r(« + è ) V2 


f\{x) 


n 2 E(-D' 


n=0 


æ 2n r(n + ^ )r( A + |) 

(2n)! r(n + À + l) 


z^ A ^ l) r " 

2/ ^ n!r(n + A + 1) V2/ 

n=i) 


2 n 


D'où l'égalité demandée. 

3.6 D'après la relation (*) f\ apparaît comme une fonction définie à droite de 0 par une série entière. 


donc 


et donc 


>. m /,(^ r ® r(A+2) 


æ-»0+ 


r(A + 1) 


Hm AMm + 1 ) = lim r(A + i)/ A (i) = l 


donc J\ ( x ) est équivalente à 


J.XXXX , „ . \ IU11 -, -, 

x— a 0+ (g) A x— a0+ r(è)r(A + è) 

1 fx\ A 


r(A + 1 ) V 2 

3.7 Une autre expression intégrale de J n , n € IN. 


2 ) x V 2 ^ V'W 2 1 

au voisinage de 0 + . 
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n+2m 


3.7.1 On a, par la formule de binôme : ( e ie - e- i9 )e~ ind = Y C k n+2m {-~l) k e {2m ~ 2k)ie , d'où : 


k = 0 


— r ( 

27T J_S 


e ie - e~ i6 )e- 


n+2m 

ine de = £ C* +2m ( 

fc =0 


l) fe r e iV™-2k)0 dQ = 

J —TT 


m 

n+2m* 


^ gll 

3.7.2 On a pour tout nombre complexe z, e z = > — -, d'où Vx > 0, Va € IR : 

^ ni 


n = 0 


(ix sin a) n ^ ( ix) n f e ia — e 


E ytx oui ex j v "> 

ni 


n = 0 


n= 0 


71 ! 


2i 


n oo 


n=0 


1 /x\ n 

rü V2. 


Y A f x ) (e îa - e~ lOL ^ n 


La fonction 


n ix sin 0—inQ 


est définie par la série de fonctions 


E;j(f )V-<r")”<r«, 

n = 0 


qui converge normalement par rapport à 9 sur [— ir, tt\, donc on peut intégrer terme à terme 


1 

2Ÿr 


i OO -, j 

—y L(*y 

2 «h k ' A2) — 


(é 9 - e- i9 )e- in9 de = 2- V ^ VI ( e ie - e - * ) " e~ inG d9 


j : 


Mais 


( e ie - e - *) ^ e - infl d0 = ]T C™(-l) m y’' e dk- 2 m-n)e dQ = | b sije- n + 2m 


m = 0 


D'où 


1 /*7T 00 

— / e* xsinO “ mO d0 = y 
2tt J „ ^ 

• /_7r m=0 


1 / X \ «+2m 


(n + 2m)! V2 


^ -j^n+2m£' 


m 

n+2m 


7» \ n JÏL / \ 2m 1 

_ \ \ / _ \ / -] \m 

2/ V 2 / m\(n- 


x 


m = 0 
oo 

E 


(ra + m) ! 

(-l) m /x^ 2m 


2 / m |p( n + m + 1) V 2 

m=0 v ' 


— Jn(ïï) 


3.7.3 D'après la relation de Chasles, on a : 


/: 


giæsin ^— inO dQ j ^ixsmO— inQ sin#— in6 


Le changement de variable t = —9 montre que 


çix sin G— inO ix sin 


D'où 
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4.1 


Ji 


âx sin O—inO 


(x) = èL 

= — r 2Re(e ixsine ~ ine )d6 

27T J 0 

i r 

= — cos(nO — x sin 9)d9. 

11 J o 


Parie IV. Application à l'étude de l'équation de Kepler 


4.1.1 La fonction tp e est de classe ^°° sur IR et Vx G IR, <p'(x) = 1 — £cos(x) > 0. De plus 

lim <p e (x) = +oo et lim <p £ (x) = — oo, donc ip e définit un ^°°-difféomorphisme de IR 

æ— >-+oo x — y — oo 

sur IR. 

4.1.2 Soit f G IR et £ G] — 1, 1[, alors il existe un seul élément de IR, noté v(t, e) tel que t = ip £ (v(t, e), 
donc on a bien : 

v(t, e) = t + e sin (v(t, e)). 


4.2 Premières propriétés de la fonction x. 

4.2.1 On a pour tout e G] — 1, 1[, 0 = 0 + e sin(0) et ir = ir + e sin(7r), donc par unicité, n(0, e) = 0 et 
v(lT, £ ) = TT. 

4.2.2 Les fonctions <p e et 1 1 — > v(t, e) sont de même parité, donc 1 1 — )• v(t, e) est impaire. 

4.2.3 Pour tout t G IR, on a v(t, e)+2n = t+2TT+£sm(v(t, e)) = t+2TT-\-£sm(v(t, £)+2 tt) = ï;(i+27r,£). 
Notons 4>(t) = v(t, e) — t, donc cj)(t + 2n) = v{t + 27r, e) — t — 2n = v(t, e) — t = <j>(t), donc la 
fonction 1 1 — y v(t, e) — t est 27r-périodique. 

4.3 Régularité de la fonction v 

4.3.1 La fonction g admet des dérivées partielles par rapport à x, £, t à tout ordre et qui sont conti- 
nues sur , donc ceci est suffisant pour dire que g est 2f°° sur l'ouvert % . 

Pour tout (x, £, t) G % , on a : 


dg_ 

dx 


(x, £, t) = 1 — £ COs(x) > 0 


4.3.2 Le théorème des fonctions implicites montre qu'on peut résoudre l'équation g(t, e, x) =0 en 
x, cette solution qu'est de classe A'°° comme g n'est autre que la fonction v de la question 4.1.2. 

4.4 Expression de x à l'aide d'une série entière de la variable e 
4.4.1 D'après théorème des fonctions implicites on a : 


et 


Donc 


d y . . 

âA' x) 


dg_ 

de 


(X,£,t) 


sinx 


dg , , 1 — £ cos x 


dx 

à*'' 1 ' 


dt 


(X,£,f) 


1 


dg . . 1 — £ COS X 


dx dx . 

dF = dt SmV - 
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4.4.2 D'après le théorème Schwarz ( v étant de classe ^°° ) et la relation de la question 4.4.1, on a : 


d_ 

de 


dv 

dt 


Slll 


d^v dv d 

d 2 v . „ dv dv . 

— sin v + n— — sm vcosv 
dtde de dt 

d 2 v dv d 

_ sm% + __ (sm%) 

d_ 

dt 


dv 

~dt 


sin n v 


4.4.3 Montrons, par récurrence sur n € IM*, que 


d n v 

de™ 


& 

dt 


\TL—1 


n— 1 


dv 

dt 


sin” v' 


Pour n = 1 on trouve le résultat démontré dans la question 4.4.1, supposons l'égalité est vraie 
à l'ordre n et montrons la à l'ordre n + 1. On a 


d n+1 v 

de n+l 


d ( d n v\ 

de \ de n J 


d ( d n ~ l 
~dê \dt n_1 


dv 

dt 


sm 


d n-l 

dt n ~ l 



dv 

dt 


sm 


d'après l'hypothèse de récurrence 


d 


\Tl— 1 


dt n ~ 

1 Vc 

d n 

f dv 

dt n 

Vdi 

d n 

/ dv 

dt n 

\dt 

d n 

( dv 

dt n 

V at 


d f dv 

— sm v 
de 


sin n v 


sin v. sin” v 


sin” +1 v 


^ d'après la question 4.4.2 


d'après la question 4.4.1 


4.4.4 Écrivons sin t sous la forme 



donc 


n+1 


(2isint)” +1 = ^ C* +1 (-l) fc e i ( n+1 “ 2fc ) t 


k = o 


et donc 


n+l 


( 2f)” +1 — ( sin” +1 t) = J2 Cn+l(-l ) k (i(n + 1 - k)r e i(n+l-2k)t 


k = 0 


et comme pour tout k G [0, n + 1], 0 < n + 1 — 2k < n + 1, alors 

n+l 


in+ 1 


A n 

-y— (sin™ +1 1) 
dt nK ’ 




k = 0 
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d'où : 


m 

-y— (sin n+1 t) 
dt n v ’ 


< (n + l) n . 


L'inégalité précédente entraîne, pour tout n G IN et tout t G IR, 


J 71 - 1 


dV 


— (sin n t) 


ni 


< 


n 


n — 1 


n! 


ieIN 


n n_1 1 

La règle de D'Alembert montre que le rayon de convergence de la série — —U' est 


donc, par comparaison, le rayon de convergence de la série en question est supérieure ou 
égale à 

e 

4.5 Expression de v faisant intervenir une série de Fourier de la variable t 

4.5.1 a n = 0 pour tout n G IN, car v £ est impaire. 

2 r 

4.5.2 Puisque v £ est impaire, b n = — / (v(t. z) — t ) sin (nt)dt. Utilisons le changement de variable 

71 J o 

t = x — e sin x, alors dt = (1 — e cos x)dx et b n devient : 

2 r 

b n (t) = — e sin x sin(nx — nesinæ)(l — ecosx)dx 

K . n 


2 r 

— / £sinx(— cosfnx — nesinx)) dx 
mr J o 


— [— esinxcosfnx — n£sinx)]n + / ecosxcosfnx — nesmx)dx 

™ \ J o 

2 r 2 r 

— / (ecosx — l)cos(nx — nesmx)dx 4 / cosfnx — n£sinx)dx 

n7r J o nn J o 

2 2 r 

— [— sinfnx — n£sinx)]n H / cosfnx — nesinx)dx 

mr mr J 0 


2 r 

— / cosfnx — ne sin x)dx. 
mr J o 


4.5.3 Le théorème de Dirichlet montre que la série de Fourier de v £ converge en tout t de IR vers 
v £ ( t ) puisque v £ est A ;1 sur IR et donc on peut écrire, pour tout t G IR : 

OO 

v(t, e) — t = b n sin (nt) 

n= 1 


et comme b n = — J n (ne ) l'égalité précédente s'écrit encore sous la forme : 


n 


Vf G IR, v(t, e) = t + 2 ^ sin(nt). 


n=l 


n 


(J. f?re )) 2 

4.5.4 Le théorème de Parseval montre que la série ^ ^ — est convergente et que 


ieIN* 


m 


E 

n = 1 


(J„(n£)) 2 1 


m 


Air 


[v £ (t)fdt = — / [x e (f)] 2 df. 
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Or ( avec le changement de variable t = x — e sin x ) 



D'où 


y- (J n (n £)) 2 

n 2 4 


4 


n=l 


Fin de corrigé 


Annexe : Une démonstration géométrique de l'équation du temps de Kepler 

Montrons d'abord qu'il existe une re- 


C S = OS— OC = ae—acosv = a(e— cosî;) 

Mais 

CS = r cos ( 7 t — 9) = —r cos(9 ) 

et donc 

(1) r cos(9) = a(cos v — e) . 

D'autre part, en vertu d'une propriété simple de l'ellipse et de son cercle concentrique. 


lation entre l'anomalie excentrique v et 
l'angle 9 sous la forme : 



a 



On a : 



donc 


rsmicM / / 

— = v 1 — E“ = => r sin(y) = a y 1 — e 2, sin v 

a sin u 

En mettant au carrée et en tenant compte de la relation (1), on trouve : 

r 2 = a 2 (e 2 + cos 2 v — 2e cos v + (1 — e 2 ) sin 2 v) = a 2 (e 2 + 1 — 2e cos v — e 2 sin 2 v 




r(l — cos 9) = a(l — e cos v) — a( cos v — e) 
o(l + e)(l — cos v) = 2a(l + e) sin 2 — 
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De même , 


= r (1 + cos 6) = a( 1 — e cos v ) + a(cos v — e) 

= a( 1 + s)(l + cos v ) = 2a(l + e) cos 2 — 

6 II + e t; 

2 V 1 -e 2 

Nous allons maintenant démontrer la vraie équation du temps de Kepler. On rappelle la deuxième loi de 
Kepler ( la loi des aires ) : 

Si S est le Soleil et M une position quelconque 
d'une planète, l’aire balayée par le segment 
[SM] entre deux positions C et D est égale 
à l'aire balayée par ce segment entre deux po- 
sitions E et F si la durée qui sépare les posi- 
tions C et D est égale la durée qui sépare les 
positions E et F. 

Voir la démonstration ci-après (Annexe2). 



2?’ 2 cos 2 I - 


D'où 


L'ellipse étudiée est la transformation affine du cercle par l'homothétie de rapport -, donc : 


Aire (S 1 AM) = -Aire (.SAP) 
a 


Or Aire(SAP) = Aire de secteur(OAP) — Aire du triangle (O SP) = -(a 2 ?; — ea 2 sinn). 


D'autre part, par la loi des aires : 


Aire(SAAf) t 


Aire de lellipse T 

où t est le temps écoulé depuis le passage au périgée ( le point A ) et T la période de la planète ( la durée 
de révolution ), donc 


Aire (SAM) = 


irabt 


D'où: 


ou encore 


2vr 

~T 


-t = v — esinn 


v = x + e sin v 


( Équation du temps de Kepler ) 


avec x = -^-t, représente la vitesse angulaire moyenne de la planète. 

Cette équation nous permet, en fonction de x ( donc en fonction du temps t ) de calculer l'anomalie 
excentrique v(t) et ensuite l'angle O(t) et le rayon r(t). 

L'équation de Kepler est transcendante, c'est à dire ne peut être résolue que numériquement. Dans le 
cas de petites excentricités, la méthode de récurrence est la plus indiquée. 

On procède par approximations successives, en posant comme première solution approchée 


Vq = X. 
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Ensuite, on substitue cette solution dans l'équation et on recommence le processus. On obtient ainsi une 
suite de solutions approchées, qui converge vers la solution exacte : 

Vo = X 

v\ = v + e sin(-uo) = x + e sin(æ) 

V 2 = x + esin(î;i) = x + esin(æ + esinæ) 


etc. 

Cet algorithme, par sa nature itérative, est particulièrement adapte a une solution par ordinateur. 

Annexe 2 : lois de KEPLER 

Théorème : Si un objet se déplace dans un champ de force central inversement proportionnel au carré de la 

distance entre le centre de l'objet et le centre de force, alors la trajectoire suivie par l'objet est une conique. 
DÉMONSTRATION : La démonstration est basée sur la deuxième loi fondamentale de la dynamique de 

NEWTON : ^ Ly :r /, = rrCf (inconnue à l'époque de Kepler). On place l'origine du système de référence 

au centre de force et on utilise les coordonnées polaires. 

On note 

M ( t ) = r(t) cos 9{t)i + r(t) sin 9(t)j 

le vecteur position de l'objet. Calculons v et 7 les vecteurs vitesse et accélération, on a : 

v(t) = [r'(t) cos 6(t) — r(t ) sin 9{t)9'(t)]i + [r'{t) sin 9(t) + r(t) cos 9{t)9'(t)\j 


et 


7 (t) = [r"(t) cos 9(t) — 2 r'{t) sin 9(t)9'{t) — r(t ) cos 9(t)(9' (t)) 2 — r(t ) sin 9{t)9"{t)]i 
+[r"(t) sin 9(t) + 2 r'{t) sin 9{t)9'(t) — r(t ) sin 9{t){9' (t)) 2 +r(t ) cos 9(t)9”(t)\j 

relation qui peut s'écrire à l'aide de vecteur radial unitaire : u = cos 9% + sin 9j et le vecteur normal 

■ ■ ^ du . . 

unitaire : n = — = — sin 9i + cos 9j sous la forme : 
dt 

^5 = Wtt) - r (^)( 0, ( i )) 2 ]^+ + r{t)9" {t))n. 


Comme la seule force en présence est radiale, alors, d'après la loi fondamental de la mécanique : 


—K 


(2) r" (t) - r(t)(9'(t)) 2 = et (3) 2 r'(t)9'(t) + r(t)9"(t) = 0 


La relation (2) s'écrit en multipliant par r : 2 rr'(t)9\t) + r 2 (t)9"{t) = 0 ou encore 
d,Q 

déduit que r 2 — = k où k est une constante ( la deuxième loi de KEPLER, 1609 ) 
dt 

La relation (2) devient donc 


dt 


,(ie_ 

dt 


= 0. On en 


d 2 r k 2 I< 

dt 2 r 3 r 2 


„ „ , , . , , 1 _ du 1 dr du du d9 

Effectuons le changement de variable u = Comme — = — - — et — = — — , et donc 

6 r dt r 2 dt dt d6 dt 

1 dr du k 

r 2 dt d9 r 2 
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ou encore 


dr du 

¥ = ~ k dë 


d 2 r 
dt 2 


d_ 

dt 

-k 


-k 


du 


de 

d 2 u k 

1¥¥ 


d 

dÔ 

e 

r 2 


du d9 

k dë\ dt 
k 2 d 2 u 
r 2 dO 2 


En remplaçant cette expression dans l'équation (4), nous obtenons l'équation différentielle linéaire d'ordre 
2 avec second membre : 

d 2 u _ K 
dë 2+u ~k 2 ' 

La solution général de telle équation s'écrit : 


u(6) 


A sin 6 + B cos 0 + 


K 

¥ 


Si nous plaçons l'axe des coordonnées polaires de telle sorte que r(0) soit minimal lorsque 0 = 0, u (e) 
devrait être maximal, ce qui donne A = 0 et, puis que u" (6) = —B < 0, donc B > 0. Ainsi 


u (e) 


b cos e + 


k_ 

¥ 


et r(9) = 


1 


r 0 k 2 k 2 

( tq = — , e = B—). Cela représente l'équation d'une conique en 

1 + e cos e K K 


K 

B COS 0+^2 

coordonnées polaires, d'excentricité e dont l'un des foyers est l'origine 


Première loi de Kepler, 1609 : Dans le référentiel héliocentrique ( le centre est le Soleil ), les centres d’inertie 
des planètes décrivent des trajectoires elliptiques dont le Soleil occupe l'un des foyers. 

En effet, d'après le théorème précédent, les trajectoires sont des coniques et comme les planètes ne 
peuvent pas s'éloigner indéfiniment du centre de Soleil, donc il y a une forte "probabilité" pour que les 
trajectoires des planètes qui gravitent autour de notre Soleil soient des ellipses. Autrement dit 0 < e < 1. 
Exemples : 


Terre e = 0, 017. Mercure e = 0, 206. Neptune e = 0, 009. 


Deuxième loi de Kepler : loi des aires, 1609 On désigne par srf Taire balayée par le vecteur position 
pendant l'intervalle de temps. 


On sait que : 



Mais comme r 2 d0 = kdt, alors 

1 f t2 k 

£/ = - kdt = ~(t 2 - h). 

2 Jt\ 2 

Le vecteur position balaie des aires 
égales en des temps égaux. 



ml26ml.tex - page 14 


Lois de Kepler : lois mathématiques décrivant le mou- 
vement des planètes autour du Soleil, formulées par 
l'astronome allemand Johannes Kepler au début du 
xviie siècle. 

Johannes Kepler établit ces lois d'après les coordon- 
nées planétaires mesurées à l'astrolabe par Tycho 
Brahé, dont il est l'assistant de 1600 à 1601. Il dé- 
couvre que les orbites planétaires ne sont pas circu- 
laires, comme le suggèrent le système de Ptolémée et 
le système de Copernic. 

Ainsi, d'après la première loi de Kepler (formulée en 
1609), les planètes gravitent autour du Soleil en sui- 
vant des trajectoires elliptiques, ce dernier occupant Première loi de KEPLER 

l'un des deux foyers de l'ellipse. 

D'après la seconde loi de Kepler (formulée également en 1609), les aires décrites par le rayon vecteur joignant la 
planète au Soleil sont égales pour des intervalles de temps égaux. 



rp2 

Selon la troisième loi de Kepler (formulée en 1619), pour toute planète gravitant autour du Soleil, le rapport Uj- 
est constant, a étant le demi-grand axe de l'ellipse correspondant à la trajectoire de l'astre autour du Soleil et T la 
période orbitale (ou période de révolution) de la planète. 

Les lois de Kepler régissent donc le mouvement de la Terre, de la Lune et des planètes, ainsi que celui des 
satellites d'un astre quelconque. C'est en s'inspirant des lois de Kepler qu'Isaac Newton découvre les lois de la 
gravitation universelle en 1687. 

Encarta 2009. 
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